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ABSTRACT 

Context. The origin of magnetism around asymptotic giant branch (AGB) stars is stih uncertain. These stars may drive an important 
dynamo, but if the magnetic energy dissipates all in X-rays, the observed X-ray luminosities are too low for a strong, dynamically 
important global field. Other explanations for the circular polarization of SiO masers in the AGB atmospheres may thus be needed. 
Aims. The interaction of the AGB wind with previously ejected matter and with planets is expected to bear complex magneto- 
hydrodynamic phenomena on a short time scale, such that strong magnetic fields can be maintained locally. Here we provide obser- 
vational evidence for the corresponding magnetic fluctuations. 

Methods. We use the circular polarization of the v = 1,7 = 2-1 SiO masers as tracer for magnetic activity. A correlation polarimeter 
allows us to simultaneously record all Stokes parameters. An SiO maser survey of 80 AGB stars was performed, out of which eight 
sources with the strongest circular polarization were selected for further monitoring. 

Results. In two AGB stars, V Cam and R Leo, we find evidence of pseudo-periodic fluctuations of the fractional circular polarization 
on a timescale of a few hours, from which we infer magnetic fluctuations of ~ 1 G. The phenomenon is rare and, if detected in an SiO 
star, restricted to a narrow range of velocities. It seems to be associated with planetary wake flows suggested by VLBI maps. 
Conclusions. While scenarios involving magnetic activity in the extended stellar atmosphere have problems to explain all observed 
features, precessing Jovian magnetospheres predict all of them without difficulty. For the case of R Leo, we constrain the orbit of the 
planet (estimated period 5.2 years), derive a stellar mass estimate of 0.7 Mq from it, and discuss the impact of planetary magnetism 
on the survival of planets. Smooth velocity variations of the fluctuating circular polarization feature are predicted as the planet moves 
along its orbit. 

Key words, masers - polarization - stars: AGB and post-AGB, atmospheres, magnetic field, planetary systems - techniques: polari- 
metric 



^ 1. Introduction the extended atmosphere of AGB stars (Nedoluha & Watson, 

_ , , , ^ , . . , . ,. 1994, Elitzur, 1996), varies in two AGB stars with a period of a 

>j Towards the end of their lifetime, low- to intermediate mass ^^^^.^ p^^^.^^^ multi-epoch observations (Diamond et al., 

^ stars undergo a phase in which they burn Helium in a shell on 2003, Pardo et al., 2004, Kang et al., 2006) of SiO masers were 

■ - top of a Carbon-Oxygen core, and Hydrogen in another shell poiarimetric, while the sampling of the polarization variabil- 

above the Helium shell. During this phase, when the stars ap- ■ ^^^^ ^^^^^ (2003) was not dense enough to detect 

pear in the Hertzsprung-Russel diagram in the upper asymp- ^^^^ .^^^^^ ^ fluctuations, 
totic giant branch (AGB), they lose, due to pulsations, at least 

half of their mass, which forms a circumstellar envelope (see Our knowledge about magnetism within the extended at- 
Herwig, 2005, for a review). The inner part of this envelope, mosphere of AGB stars relies on circulai- polarization measure- 
also called extended atmosphere, is expected to bear complex ^^nts of SiO masers at 1.5 to 7 AU distance from a star with a 
magneto-hydrodynamic phenomena, due to the interaction of the typical radius of 0.7 to 1 AU. Maser theory implies molecular ex- 
wind with the previously ejected matter and with planets (Struck "Nation in dense pockets of gas and amplification of light in nar- 
et al., 2004, 2002, Struck-Marcefl, 1988). As in the solar sys- ^w tubes. Assuming that the circular polarization is caused by 
tem, where space weather changes on timescales of hours (e.g. *e Zeeman eff-ect in the non-paramagnetic SiO molecules in the 
Prange et al., 2004), fluctuations of the magnetic field about ^^^es, the reported fractional circular polarizations of up to 9 % 
a mean value can be expected due to this interaction, but the for the y = 1 , 7 = 1 - transition (y and J are the vibrational and 
observational evidence is still lacking. Here we wiU show that rotational quantum numbers, respectively) yield line-of-sight av- 
for a narrow range of velocities the circular polarization of SiO ^""^8^^ magnetic flux densities of up to 100 G (Bai-vainis et al., 
masers, generaUy accepted as a tracer of the magnetic field in 1987). More recent observations (Herpin et al., 2006) of the 
V - 1,7 = 2-1 masers (excited in gas layers very close to 



Send offprint requests to: H. Wiesemeyer those where v - 1,7= 1 - is excited) also show the presence 

* On leave to Instituto de Radioastronomi'a Milimetrica, Granada, of circular and Unear polarizations in many stars. These observa- 
Spain. tions result in lower magnetic flux density estimates in the range 
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1 to 15 G, again assuming the circular polarization is due to the 
Zeeman effect. In the outer envelope, H2O and OH masers re- 
side at ~ 100 to 400 AU, respectively ~ 1000 to ~ 10000 AU, 
from the star. The Zeeman effect of H2O masers has been used 
(Vlemmings et al., 2002) to estimate, by extrapolation assuming 
a solar-type field topology, a magnetic flux density at the stellar 
surface of 100 G. In summary, it is now observationally evident 
that AGB stars maintain a magnetic flux throughout their enve- 
lope. Blackman et al. (2001) showed that AGB stars can generate 
a magnetic field via dynamo action at the interface between the 
rapidly rotating core and the extended convection zone. The im- 
portance of the magnetic field for global envelope dynamics is 
still a matter of debate, though (e.g. Soker, 2006). 

Once the SiO maser spots are formed (often distributed along 
arcs around the star, hereafter referred to as maser shell, see e.g. 
Cotton et al., 2008), they are subject to magnetospheric events 
which are known to be rapidly variable if we refer to the so- 
lar system. Here we provide first evidence for such fluctuations 
in the atmospheres of two AGB stars by frequent polarization 
sampling. The polarization monitoring by Glenn et al. (2003) 
was only sensitive to slow variations and therefore to a long- 
term readjustment of the magnetic field (on a timescale of sev- 
eral months). 



2. Observations and data analysis 

2.1. Instrumentation and observations 

The bulk of our observations was done in May 2006 with the 
IRAM 30m telescope at Pico Veleta, which is equipped with 
dual-polarization receivers. To confirm the rapid polarization 
variability discussed here, Mira-type star RLeo was observed 
during additional 4 hours on 2008 August 01. The receiver pair 
operating at 3mm wavelength was tuned to the v = 1,7 = 2-1 
transition of SiO at 86.243 GHz and made coherent by us- 
ing the same local oscillator reference. The signals were an- 
alyzed with XPOL (Thum et al., 2008), a correlation spec- 
trometer enhanced with cross-correlation products for the sig- 
nals from the orthogonally polarized receivers, at 39.0625 kHz 
or 0.136kms ' channel spacing within a 16.2MHz bandwidth. 
Receiver temperatures varied between 36 K and 64 K, and sys- 
tem temperatures between 70 K and 235 K, resulting in an an- 
tenna temperature (T^) noise of cr„„s = 36 - 100 mK across the 
spectral baseline. The conversion factor from T* to the flux den- 
sity scale is 6 Jy/K. The phase difference between the receivers 
was measured with a polarizing grid mounted in front of the cold 
load (at about 70 - 80 K effective temperature) of the calibra- 
tion unit. Whenever the temperature scale was calibrated (ev- 
ery 10 to 20 minutes), the signal from the polarizer, with well 
known properties, was also observed and compared to the unpo- 
larized calibration load at ambient temperature. The phase cor- 
rection was then applied in order to attribute the measured real 
and imaginary part of the cross correlation product to the Stokes 
U parameter (in the reference frame of the telescope's Nasmyth 
focus) and to the Stokes V parameter (positive for a right hand 
circularly polarized signal, according to lAU convention). Small 
residual phase errors result in a leakage of the Stokes U signal 
into Stokes V, a subtle effect only worth worrying about when 
Stokes U is strong (which is the case here). The scaled (~ 3 %) 
copy of the Stokes U signal thus contaminating the measured cir- 
cular polarization can be removed, since it is subject to the par- 
allactic rotation of the polarization vector, which easily shows in 
the data. The remaining Stokes V signal is intrinsic, the contri- 



bution of the telescope and the receiver cabin optics being neg- 
ligible (-0.03 + 0.12 % on the optical axis). 

77 AGB stars were observed, out of which 62 are Mira stars, 
10 are semiregular variable stars, and 2 red supergiants. Eight 
objects with the strongest circular polarization were retained 
for dense monitoring (twice per hour). They are all Mira stars, 
in agreement with the finding of Herpin et al. (2006) that this 
source class tends to have the strongest linear and circular polar- 
izations, and is therefore well suited for our aims. 

2.2. Data analysis 

The temperature and phase calibration was done with the MIRA 
raw data reduction software'. The subsequent processing of the 
Stokes spectra, described in the following, was done with the 
CLASS and GREG softwar^B- Due to the irregular sampling, a 
determination of the spectral power density (SPD) of the time 
series using Fourier transform techniques would yield poor re- 
sults, and preference was given to the SPD estimate provided by 
the Lomb technique (Lomb, 1976, see also Press et al., 1994). 
This method only uses the measured data without any prior inter- 
polation to a regular sampling function, and does a least-square 
fitting to the harmonic contents of the time series. A detection 
of quasi-periodic fluctuations of the fractional circular polariza- 
tion (hereafter pc) needs to pass three critical tests. First, we 
do not expect a significant fluctuation in the Stokes U residuals 
left after the subtraction of the stationary (within the timescales 
considered here) linear polarization feature. Second, the oscilla- 
tion of Pc was tested against the null hypothesis that the data are 
not periodic, but random noise with a Gaussian distribution; this 
test yields a false-alarm probability. Third, because a restricted 
number of random samples cannot be strictly Gaussian, Monte- 
Carlo simulations, with random samples contemporaneous with 
the observed features, were used to confirm the significance of 
the false-alarm probabilities. 



3. Results 

Only two stars, the Mira variables V Cam and R Leo (Fig. 1) 
passed the three tests and show clear evidence for a quasi- 
periodic modulation of their fractional circular polarization pc 
(Fig. 2), of peak-to-peak amplitude 2.8 % and 3.7 % for V Cam 
and RLeo, respectively (Tab. 1). The observations were made 
~ 40 days after the optical maximum of V Cam (pulsation 
period 522 d) and ~ 30 days after that of R Leo (period 
310d), in an expansion phase of the SiO maser shell. Within 
a ~ 0.4km s"' narrow range of radial velocities centred at 
7.5 km s"' and 4.4 km s"' , we detect periods of 5*^4 + d?l and 
6*^3 + 0'?3, for V Cam and R Leo, respectively (Tab. 1). The phe- 
nomenon is both rare and localized in velocity space (and thus 
in circumstellar space, due to the velocity structure). The likeli- 
hood of misinterpretation of a time series of Gaussian noise is 
< 4 % and < 16 % for V Cam and R Leo, respectively (Fig. 3, 
left panels). We note that the Lomb periodograms of the Stokes 
U residuals do not show any significant periodic signal (Fig. 3, 
right panels). 

The Monte Carlo tests confirm the significance of the results. 
For 10^ runs, only 0.28 % of the Lomb periodograms mimic a 
periodic signal with < 4 % false-alarm probability (V Cam), re- 
spectively 2.9 % a periodic signal with < 16 % false-alarm prob- 
ability (R Leo). We therefore conclude that it is unlikely (for 
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(1) 

(2) 
(3) 
(4) 
(5) 





V ^^am 


K i^eo 




u.uo 


0.10 


Oscillation period of circular polarization'^^' 


(5.4±0.1)h 


(6.3 ± 0.3) h 


False-alarm probability'^' 


4% 


16% 


Radial velocity of maser spots'"*' 


7.5kms-' 


4.4 km s"' 


Mean circular polarization 


-0.3 % 


-1-2.0% 


Mean magnetic flux density in maser'^' 


~ 100 mG 


~ IG 


Peak-to-peak variation of fractional polarization 


2.8% = 5.6a-,„s 




Fluctuation of magnetic field*'' 


~ IG 


~ IG 



With respect to the maximum of the optical lightcurvc, 
for V Cam the second harmonic of the modelled rotation period, 
with respect to the null hypothesis (data are Gaussian noise), 
with respect to the local standard of rest, 

order-of-magnitude estimate based on the Zeeman hypothesis, and for a linewidth of 1 km s"' . This corresponds to a maser at 3 rj and a 
planetary magnetic dipole field eight times stronger than Jupiter's (Wiesemeyer, 2008). 





Fig. 1. Spectra of total power (Stokes I, i.e. polarized plus unpolarized flux density) and of fractional circular polarization. Grey- 
shaded histogram: Stokes I flux density in Jansky vs. vis, (the radial velocity with respect to the local standard of rest). Dots with 
errorbars: fractional circular polarization pc,in% with respect to Stokes I. Left: Mira-type star V Cam. Right: Mira-type star R Leo. 



V Cam extremely unlikely) that random noise creates the ob- 
served fluctuations. If they were produced by an oscillation of 
the telescope tracking, moving the polarized sidelobes across 
the line-of-sight towards the source, all velocity channels would 
show fluctuations with the same period, which clearly is not ob- 
served and can therefore be safely excluded. 



Variations of the circular polarization of 4 % peak-to-peak 

amplitude were found again in August 2008 (Fig. 4), but now 
within two and wider velocity intervals, from visr = -3.6 to 
— 2.7kms~^ and from -0.5 to 0.7kms~' . A meaningful Lomb 
analysis as in the case of the 2006 data was not possible for these 
data, because the observing interval was too short to measure a 
full fluctuation period. The data (Fig. 4 left and middle) suggest 
a period of about 2 hours with a guessed uncertainty of about 
0.5 hours. It is interesting to note that in channels of relatively 
strong Stokes I emission, no significant polarization fluctuation 
was found (Fig. 4, right), while now also the linear polarization 
seems to fluctuate at about visr= 0.6 km s"' , where Stokes I is 
strong. Again, this rules out artifacts as the reason for the fluc- 



tuations. If they were of instrumental origin, the whole spectral 
band would be affected. 



4. Discussion 

The fractional circular polarization is a more reliable magnetic 
tracer than hnear polarization, especially for spatially unresolved 
observations: our telescope beam averages the polarization sig- 
nal from all maser spots at the same radial velocity, hence the 
orientation of the linear polarization is either parallel or per- 
pendicular to the maser shell (e.g. Diamond & KembaU, 2003, 
Desmurs et al., 2000). This tends to decrease the beam-averaged 
hnear polarization, except if it is dominated by few, very strong 
maser spots close together (which may explain the linear polar- 
ization fluctuations showing up in the August 2008 observations 
of R Leo, Fig. 4, middle). Furthermore, linear polarization can be 
entirely produced without magnetic fields by anisottopic pump- 
ing (Western & Watson, 1983, Asensio Ramos et al., 2005), 
though polarization angle swings by 90° in the afore mentioned 
resolved observations are explained by changes of the magnetic 
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VCAM 

Channel 163, V|,r [km/s] = 7.5 



RLEO 

Channel 201, v,,, [knn/s] = 4.4 




20 25 
UTC [hours], 22 May 2006 




UTC [hours], 22 May 2006 



Fig. 2. Time series of polarization measurements. From top to bottom: Stokes I in Jansky, the fractional polarizations pl (linear) 
and pc (circular), and the polarization angle (i.e. the position angle of linear polarization, in degree E from N), for the SiO maser 
spots from V Cam at vis,-= 7.5 km s"' (left) and R Leo at visr = 4.4 km s"' (right). The continuous heavy lines in the plots for pc 
show model results (Wiesemeyer, 2008) for a saturated maser in the equatorial plane, at 3 rj from the planet, with a magnetic dipole 
field of eight times Jupiter's and misahgned with the rotation axis (in the sky plane) by 10°, and rotation periods of 10'?8 (V Cam, 
left) and ^3 (R Leo, right). 



field direction with respect to the propagation direction of radi- 
ation in the maser (Goldreich et al., 1973). By contrast, circular 
polarization is enhanced in magnetized maser spots, either di- 
rectiy via the Zeeman splitting (e.g. Watson & Wyld, 2001), or 
indirectly via birefringent conversion of linear into circular po- 
larization due to changes of the magnetic field topology along 
the maser slab (Wiebe & Watson, 1998, for unsaturated masers), 
or due to maser saturation, rotating the quantization axis from 
the direction of the magnetic field to that of the radiation prop- 
agation in the maser (Nedoluha & Watson, 1994). In the latter 
case, pump conditions would be required not only to vary in time 
but also quasi-periodically, an event which we cannot exclude 
but seems unlikely on short time scales. In turn, this would pro- 
duce a varying Stokes I flux density, which is not observed. We 
therefore suggest that the intraday fluctuations of pc are of mag- 
netic origin. They are quasi-periodical, not stochastic, and rare 
and only at one velocity, which means that an ordered magnetic 
field structure with smooth gradients in time and space is locally 
maintained within the dense AGB winds of R Leo and V Cam. 

4.1. A case for precessing Jovian magnetic fields ? 

The magnetopheres of Jovian planets have already been pro- 
posed to explain several features of SiO maser polarization 
(Struck et al, 2002, 2004), thus providing a local explanation 
for the strong magnetism in the atmosphere of AGB stars. If the 
magnetic fields in which SiO masers are embedded, of typically 
10 G, were global, one would assume, at the stellar surface, a 



magnetic field of about 100 to lOOOG, depending on assump- 
tions on the gradient of the magnetic flux density. Such a field 
could be generated by a dynamo between the rotating core and 
the stellar convection zone (Blackman et al., 2001). However, 
we now have to explain (see item 3 below) why Mira stars were 
not detected (Haisch et al, 1991) by the Rosat all-sky soft X-ray 
survey, which is sensitivity-limited to X-ray luminosities above 
typically 2.4 x 10^^ erg s"' . As a matter of fact, Kastner & Soker 
(2004a) find that the X-ray luminosity of the single Mira stars 
TXCam and TCas is below the threshold expected for a dy- 
namically important magnetic field. Their conclusion relies on 
the assumption that the X-ray luminosity equals the kinetic lu- 
minosity of the stellar wind powered by magnetic fields. Since 
there may be ways to dissipate magnetic energy other than by 
X-ray emission from accelerated particles, this is still an open 
issue. 

We now summarize the arguments supporting the plausible 
planetary origin of our detections. 

1. The rareness of the phenomenon can be explained by the 
number of Jovian planets around AGB stars, and by the re- 
quirement that the gain region of the maser needs to cross 
the magnetosphere where it is dominated by the dipole com- 
ponent. At present we cannot distinguish which explanation 
prevails - the possible rareness of Jovian planets in AGB 
atmospheres due to infall of planets into the star (Villaver 
& Livio, 2007, Nordhaus & Blackman, 2006), or the un- 
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V|g^ = 7.5 km/s 



0.2 0.4 0.6 

frequency [1/hours] 




0.2 0.4 

frequency [1/hours] 



R Leo 



V|3^ = 4.4 km/s 
0.1 




R Leo 



v,.. = 4.4 km/ s 




0.5 

frequency [1/hours] 



Fig. 3. Lomb periodograms of the polarization time series (estimates of the spectral power density by least-square fits of their 
harmonic contents, Lomb, 1976), normalized by the variance of the samples. Top left: Lomb periodogram of pc from V Cam at 
visr - 7.5 km s"' . Top right: same for the residuals of Stokes U (normalized by Stokes I, after subtraction of the stationary linear 
polarization). Bottom left and right: as top panels, but for R Leo at visr = 4.4 km s ' . The horizontal dashed lines indicate the spectral 
power density at which a periodic signal can be mimicked by Gaussian noise with a probability given by the labels ("false alarm 
probability"). For a given source, the plot scale for the Lomb periodograms of pc and the Stokes U residuals is the same, to make 
the comparison of the significance of the peaks easier. 



favourable location of the SiO maser with respect to the pre- 
cessing magnetosphere. 

2. Both in V Cam and R Leo, the signature is statistically sig- 
nificant only at one velocity. We interprete this velocity as 
the current line-of-sight velocity of a planet on its Keplerian 
orbit. 

3. Dedicated models (Wiesemeyer, 2008) of saturated SiO 
maser emission in a precessing magnetic dipole field show 
that a dipole moment eight times stronger than that of 
Jupiter's magnetic field is sufficient to provide the ~ 1 G 
magnetic flux density at 3 rj distance, required to explain the 
peak-to-peak variation of pc (Tab. 1). This reconciles maser 
polarization and the lack of an X-ray luminosity expected for 
a stellar corona. Saturated maser action is not required here, 
if the weak circular polarization of unsaturated masers is en- 
hanced by the combined efl'ect of anisotropic pumping and 
magnetorotation (Wiesemeyer, 2008). 



4. The fluctuations are quasi-periodic, not stochastic. The fluc- 
tuation periods derived from the Lomb periodograms of the 
Pc fluctuations in R Leo and V Cam can be modelled as the 
first and second harmonic, respectively, of the planetary rota- 
tion period which is &.?> and \&.'S>, respectively (cf. Jupiter's 
rotation period of 9'?9, and Saturn's of 10'?5). A misalign- 
ment of the dipole axis with the planetary rotation axis (10° 
in the models, 9?6 for Jupiter) naturally explains both the am- 
plitude and the period of the fluctuation, due to the ordered 
structure of the magnetic field. This implies that the mag- 
netosphere is protected by a Chapman-Ferraro type magne- 
topause, where ram pressure equals magnetic pressure, at 
about lOrj (Jupiter radii) substellar distance (according to 
Struck-Marcell, 1988). To date it is still uncertain why our 
RLeo data from 2008 yield a different period. This may be 
partly due to the incompleteness of the data, partly due to the 
proper motion of the maser spots, which may cross 10 rj in 
10 hours and thus modify the pattern expected for a station- 
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Fig. 4. Same as Fig 2, but for R Leo at three selected velocities (observations from 2008 August 01). Left: Vkr = -2.8 km s ' . Middle 

: visr=0.6kms ' .Right: visr= -5.3 km s"' 



ary fluctuation. Any intrinsic origin (changes of the angular 
momentum of accretion onto the planet, or of the spin of 
the latter) would be most difficult to model. Likewise, it is 
still unclear why in 2008 polarization fluctuations appeared 
in two well separated velocity ranges. 

5. If the pseudo-periods of the fluctuations were only due to 
the proper motion (~ lOkms ' ) of SiO maser spots (with 
~ 0.1 AU diameter) across a quasi-stationary magnetic per- 
turbation, this would result in periods which are two orders 
of magnitude longer than the observed ones. In turn, this 
means that the magnetic fleld diffuses through the gas of the 
maser slab, i.e. the ambipolar diffusion speed exceeds the 
maser velocity. A Jovian magnetic field both provides the re- 
quired strength and curvature radius of the field lines (Struck 
et al., 2002), unlike a magnetic field of stellar origin. 

6. Spatially resolved observations (Cotton et al., 2008) of the 
V = 1,7 = 1-0 SiO masers in R Leo show, at ~ 4 - 
7.9 km s"' radial velocity (with the largest redshifted veloc- 
ity closest to the star), a remarkable elongated feature sug- 
gestive of a planetary wake flow as modelled by Struck et al. 
(2004). The corresponding feature in the v = 2, / = 1-0 
transition is even double-lobed. Similar features were found 
by Cotton et al. (2004) and Soria-Ruiz et al. (2007) at various 
circumsteUar positions and difl'erent epochs. 

4.2. Kinematical and astrometric evidence for a planet 
orbiting around R Leo 

We deliberately consider here the presence of elongated fea- 
tures (see item 6 above) as the possible signature of a plane- 
tary wake although explanations involving the theory of maser 
formation cannot be ruled out. In the frame of this working hy- 
pothesis we now will use the available kinematic and astromet- 
ric information to illustrate the potential interest of combining 
SiO polarization monitoring observations with contemporane- 
ous VLBI imaging of SiO masers. The VLBA observations of 



Cotton et al. (2008) were in September 2004, 638 days after sim- 
ilar observations by Soria-Ruiz et al. (2007) in December 2002, 
while we observed in May 2006, 624 days after Cotton et al. 
(2008). Remarkably, the v=l,/=l-0 map of Soria-Ruiz 
et al. shows, at about visr= -7kms"' (with respect to the local 
standard of rest, hereafter LSR), a similar elongated maser fea- 
ture pointing radiaUy away from the star towards the South-East, 
and another one towards North-East, with velocities peaking at 
visr = 2.5 kms"' and -2kms"' . This suggestion is confirmed by 
inspection of the VLBA map taken in August 2001 by Cotton 
et al. (2004), again showing a quasi-radial aUgnment of maser 
spots, now towards the north. The repeated occurence of such 
features, systematically rotated counter-clockwise from North, 
strengthens the suggestion of a wake flow associated with the 
same planet. However, the combined use ofv = 1,7 = 2-1 
polarization monitoring with the v = l, 7=1-0 VLBI astrom- 
etry raises new questions. Soria-Ruiz et al. (2007) show that in 
R Leo the former maser spots do not coincide with those corre- 
sponding to the latter transition. The characteristic ring radii of 
the respective maser shells are 29.2mas and 33.8mas, for the 
7=1-0 and 7 = 2-1 transitions, whereas coUisional and 
radiative pumping models (see Soria-Ruiz et al. for further ref- 
erences) suggest that for our purposes the maser spots of these 
transitions should be sufliciently close to each other. With this 
inconsistency in mind a lower Umit to the astrometric error in 
the comparison of our polarization data with the VLBA posi- 
tions can therefore be estimated from the difference between the 
ring radii, 4.6 mas or 0.5 AU (for 113.5 pc distance, Fedele et al., 
2007, further references therein). The actual astrometric error 
may be higher, because the absolute astrometry of these VLBI 
data from different epochs is far less well constrained than the 
relative astrometry for a single observing date. Despite this limi- 
tation we use here the available combined kinematics and astro- 
metric information to constrain the planetary orbit. The velocity 
in the wake ffows is usually decelerated, the maximal blue- re- 
spectively redhifted velocity should therefore be closest to the 
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planet. A sine fit (Fig. 5, top) yields a period of 5.2 years and 
a stellar velocity of visr = 1 .7 km s~' , but more observations are 
needed in order to assess the significance of the fit. Our 2008 
data are best described by the velocity of the more blueshifted 
polarization fluctuation. The fitted stellar velocity is between the 
value derived by Bujarrabal et al. (1989) from thermal molec- 
ular lines in the circumstellar envelope (-0.5 km s^' ) and the 
catalogue value of 13.4km s-^ (GCRV, Wilson, 1953). Using the 
method discribed in Appendix A, this velocity fit allows us to 
determine the radius of the planetary orbit to 23.8 mas which 
corresponds to 2.7 AU, at a distance of 1 13.5 pc. The inclination 
of the orbital plane with respect to the plane of the sky is i = 34°. 
The result of the least-square fit to the available astrometric in- 
formation is shown in Fig. 5 (bottom). The linear polarization 
measured by us is predominantly parallel to the radius vector 
calculated from the model fit. Due to uncertainties in our method 
and in the astrometric positions of the VLBI spot maps (no ab- 
solute uncertainties are derived in the works of Cotton et al. and 
Soria-Ruiz et al.), it is not too surprising to observe a discrep- 
ancy between the observed and modelled positions, especially 
for the 2002 Dec 07 VLBA observation (Soria-Ruiz, 2007). 

As for the August 2008 data in RLeo which show polar- 
ization fluctuations at two different velocities, contemporaneous 
VLBI observations would have been required to perhaps fur- 
ther suggest that the second velocity interval could be due to 
a wake flow and therefore to another planet. However, it seems 
too speculative at this stage to suggest that each elongated VLBI 
maser feature could be associated with a planet while jet- like 
features could as well be interpreted as ejections of matter (see 
e.g. Cotton et al., 2008). Our fit leads to a stellar mass estimate 
of 0.7 Mq, below, but not unreasonably far from the value in- 
ferred from measurements of stellar parameters by infrared K- 
band interferometry (Fedele et al., 2005), namely 1.0 - 1.2 Mq, 
for an equivalent non-pulsating star of the same luminosity. The 
discrepancy can either be due to observational error or model un- 
certainties. As a matter of fact, the R Leo maps of Cotton et al. 
(2008) mentioned above show several elongated maser features, 
and it is unlikely that each is associated with a planet (Cotton et 
al. rather suggest jet-hke ejections of matter, see section 4.3.6 be- 
low). The conclusion that the polarization fluctuations observed 
at two velocities in R Leo in 2008 are due to the presence of two 
planets would therefore be speculative, even if contemporaneous 
VLBI data were at hand. 



4.3. Are alternative scenarios viable ? 

Boldly extrapolating the rich magnetic phenomena in the solar 
system to late-type stellar systems, alternative scenarios have to 
be seriously examined. One might suggest magnetic activity due 
to stellar companions, coronal loops, magnetic clouds or Parker 
instabilities in shocks forming in response to steUar pulsations. 



4.3. L Binarity 

Here the mass transfer from the AGB star to a close companion 
generates flares and related magnetic fluctuations, which are not 
necessarily located in the atmosphere of the former (e.g. the Mira 
AB system, Kastner & Soker, 2004b). However, neither V Cam 
nor R Leo are known for having stellar companions, otherwise 
their X-ray luminosities would be higher. As a matter of fact, 
Tatebe et al. (2006, IR interferometry at 1 1.5 fim) do not see any 
obvious asymmetry in the envelope of RLeo that could hint at a 
companion (see also Gatewood, 1992). 
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Fig. 5. Top: Keplerian velocity fit of maser features associated 
with the suggested planet (LSR velocities at the tip of the plane- 
tary wake flow, and velocities of polarization fluctuations). Time 
is in years since 2001 Aug. 04 (observations of Cotton et al., 
2004). Other references: Soria-Ruiz et al., 2007 (2002 Dec. 07), 
Cotton et al., 2008 (2004 Sep. 05), this work (2006 May 22 and 
2008 Aug. 01). Bottom: available astrometric information and 
least- square fit of the derived planetary orbit (Appendix A), for 
d = 113.5 pc stellar distance. VLBI positions are indicated by 
black dots. The planetary orbit is indicated by the black line, 
crosses mark the fitted position at the VLBI epochs. Empty cir- 
cles give the positions of the planet at the epochs of our polariza- 
tion measurements. The linear polarizations measured by us and 
their position angles are shown as bars (the bar in the top left 
corner indicates a fractional linear polarization of 30 %). The 
grey-shaded disk indicates the measured A 2.2 fim uniform-disk 
diameter at stellar phase 0.08 (in April 2001, Fedele et al., 2005). 
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The case of the 2003 soft X-ray outburst in the Mira AB sys- 
tem is, interestingly, not related to the binarity of the system. The 
spatially resolved Chandra observations (Karovska et al., 2005) 
clearly show that it occured within the envelope of Mira A, and 
it would be worthwhile to estimate whether it was triggered by 
burst-like accretion onto a Jovian planet, as suggested by Struck 
et al. (2004). As a matter of fact, the position of the X-ray flare 
may coincide with a linear SiO feature seen by Cotton et al., 
(2004). 



4.3.2. Coronal loops 

Coronal loops are generated by the difi'erential rotation in the 
sub-photosphere, together with the dynamo (Blackman et al., 
2005) at the interface wifli the core. If the convection zone 
extends to at least 80 % of the stellar radius (Holzwarth & 
Schiissler, 2001), magnetic flux loops are trapped there, because 
magnetic stresses at the tips of the loops act against convec- 
tive buoyancy. This also explains the lack of a coronal X-ray 
luminosity, and why the magnetic fluctuations are rare. One 
may suggest coronal loops to occasionally extend to the site of 
the SiO masers. If there were coronal flux loops, the timescale 
for magnetic fluctuations would be too long, though. This can 
be demonstrated in the following way: Assuming (by analogy 
with the sun) that the size of a loop corresponds to that of a 
super-granulation cell, this yields 0.2 AU, for 1 AU stellar ra- 
dius and 400 convection cells in the stellar sub-photosphere 
(Schwarzschild, 1975), in contrast to 2 X 10^ for the sun. 
Interferometric K-band size measurements (Fedele et al., 2005) 
yield a Rosseland radius of 1.6AU at phase 0.08, close to the 
brightness maximum in April 2001 (cf. our measurements at 
phase 0.08 after the brightness maximum in 2006). Our orbital 
radius estimate of 2.7 AU thus infers that the maser is 1.1 AU 
above the stellar photosphere. Assuming a photospheric rotation 
speed of ~ 1 km s"' (a higher rotation speed is only possible for 
AGE stars with at least a brown dwarf companion, Soker 2006), 
a corotating flux loop would need ~ 100 days to cross a maser 
spot of 0.1 AU diameter. This is two orders of magnitude above 
our fluctuation periods, which makes it highly unlikely that coro- 
nal loops can generate the observed hour-scale magnetic fluctu- 
ations. The scenario involving trapped magnetic flux loops also 
demonstrates that from the weak X-ray luminosities one cannot 
infer the absence of an efficient stellar dynamo - the magneti- 
cally accelerated particles simply remain within the convective 
envelope and do therefore not show in X-ray observations. 



4.3.3. Magnetic clouds 

As coronal loops, magnetic clouds originate from the magnetic 
activity on the surface of the stellar photosphere. If coronal loops 
are absent, due to a lack of the latter, this also holds for magnetic 
clouds of the type encountered in the solar wind (e.g. Bothmer & 
Schwenn, 1998). Should magnetic clouds traverse the extended 
stellar atmosphere, we can provide a rough estimate of the time 
scale associated with the corresponding magnetic fluctuations. 
The size of magnetic clouds in the solar system is comparable 
to that of SiO masers (0.1 AU), and their magnetic field vector 
also rotates on a timescale of ~ 10 hours. However, solar sys- 
tem magnetic clouds expand with the speed (~ 400 km s"' ) of 
the magnetically driven solar wind. Since the latter is an order 
of magnitude above that of an AGB wind and SiO maser proper 
motions, the timescale does not fit - the resulting polarization 



would vary too slowly (especially if the magnetic cloud is decel- 
erated by the dense ambient gas). 

4.3.4. Parker instabilities 

In order to produce polarized maser features of the observed life- 
time (months to years, Glenn et al., 2003), the latter needs to 
be shorter than the ambipolar diffusion timescale, i.e. the mag- 
netic flux has to be coupled to the masing matter. The fragments 
forming as a consequence of Parker instabilities (Hartquist & 
Dyson, 1997) fulfil this requirement under the typical condi- 
tions in the extended atmosphere of AGB stars. Since Parker in- 
stabilities arise from a local indentation orthogonal to the field 
lines, some short-term variations of the magnetic field compo- 
nent along the line of sight may be expected if the ambipo- 
lar diffusion timescale of these fluctuations exceeds that of the 
maser lifetime and proper motion. The resulting variability of 
PC should globally affect all maser spots of an SiO star if they 
commonly arise from fragments formed by this mechanism, in 
contradiction with our observations. Furthermore, the mass of 
fragments forming in response to Parker instabilities scales with 
magnetic pressure, and the strongest fluctuations would be ex- 
pected from the strongest masers, which has not been observed 
(cf. Figs. 2 and 4). 

4.3.5. Alfven waves 

Our discussion would not be complete without considering 
Alfven waves (Alfven, 1942, see e.g. Sturrock, 1994) as a pos- 
sible origin of rapid magnetic fluctuations. However, our knowl- 
edge of hydromagnetic wave propagation in the dense AGB at- 
mospheres is poor. The relevant time scale here would be 

tA = D/VA (1) 

where D is the diameter of a maser slab, and va is the Alfven 
speed for the propagation of a hydromagnetic wave, given by 
(cgs units) 

B 

VA = (2) 

where in the weak-coupling regime p is the ion mass density, 
and in the strong coupling regime (i.e. when collisions between 
ions and neutrals are important) the total mass density of neutral 
and ionized particles (Zweibel, 1989). Which regime holds in the 
atmosphere of AGB stars is uncertain. If the ambipolar diffusion 
speed is close to the velocity of the AGB wind, as estimated by 
Soker (2006) and Hartquist & Dyson (1997), strong coupling 
holds, and we estimate the Alfven speed to 

lOG 1 

Va = ~ 200 kms"^ . (3) 

V47r2.10-i4gcm-3 

for a 10 G magnetic field and a hydrogen number density of 
lO'^cm"'. Thus, the timescale for an Alfven wave crossing a 
maser slab of 0. 1 AU diameter is about 20 hours, not unrea- 
sonably far from the timescales given in Table 1. However, we 
would like to stress that the parameters used for the estimate of 
both the ambipolar diffusion and Alfven speeds are much uncer- 
tain. 

4.3.6. A precessing jet ? 

The role of a planetary or brown-dwarf companion in the shap- 
ing of circumstellar envelopes has been modelled by Nordhaus 
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& Blackman (2006). Since in their scenarios the companion is 
already engulfed in the convective envelope of the star, it cannot 
be detected anymore by polarimetry. The deposit of the plan- 
etary angular momentum leads to a spin-up of the star and by 
consequence to the ejection of matter by equatorial mass loss 
or, if assisted by the stellar magnetic field, to a bipolar jet. Is 
it conceivable that we observed rather the polarization signature 
due to the precession of such a jet ? If the precession is caused 
by the misalignment between the jet axis and the stellar rotation 
axis, we can estimate the timescale of the phenomenon as in the 
case of coronal loops (section 4.3.2.) except that we have to ac- 
count for the spin-up of the star Nordhaus & Blackman showed 
that an infalling brown dwarf of 0.05 Mq may increase the stel- 
lar rotation rate by an order of magnitude. The fluctuation of the 
circular maser polarization due to a crossing jet would therefore 
occur on a timescale of at least 10 days if the cross-sectional area 
of the jet perpendicular to its axis is smaller than the size of the 
maser spot, or longer otherwise. This is still one to two orders of 
magnitude above the fluctuations observed by us. Furthermore, 
to leave a signature only in a few SiO masers the opening angle 
of the jet cannot be too large. 

5. Conclusions 

In summary, in our opinion none of the non-planetary alterna- 
tives satisfactorily explains the present observations of RLeo 
and V Cam, while a single scenario, namely precessing Jovian- 
type magnetospheres, are consistent with all of them. Since the 
discovery of a planet around a star in the post-red giant phase 
(Silvotti et al., 2007), we know that planets can survive in such 
an environment. Whether a given planet survives the AGE phase 
of its parent star depends on several factors (Villaver & Livio, 
2007). Evolution of the orbit is the main one - depending on its 
distance from the star, the planet is either dragged inwards by 
tidal friction and evaporated, or repelled due to stellar mass loss, 
while for Jupiter-sized planets the gain of mass due to Roche 
flow from the star or accretion from the stellar wind barely affect 
their orbits. If the suggested planet around R Leo, one of the best 
studied Mira stars, of solar mass and at 11 3.5 pc distance (Fedele 
et al., 2005), does not spiral in and evaporate during the remain- 
der of the upper AGB phase, it will survive the subsequent plan- 
etary nebula stage against evaporation (Villaver & Livio, 2007) 
only if it has at least two Jupiter masses and orbits at least at 
~ 3 AU. We note that the fluctuating SiO maser feature in R Leo 
has an orbital radius of 2.7 AU, close to the critical distance. We 
have demonstrated that the orbital elements of its host planet can 
be determined by measuring the velocity of the fluctuation and 
the position of the corresponding maser spot. The orbital period 
of the planet is 5.2 years, which provides an estimate of the stel- 
lar mass of 0.7 Mq. This agrees with mass estimates found in 
the literature, of about 1 Mq for R Leo. On the long term, our 
method has the potential to reveal, together with an independent 
estimate of the stellar mass, sub-Keplerian orbital motion and 
thus migration of the planet through the atmosphere towards the 
star. This requires correspondingly accurate observations. The 
situation may be complicated by the presence of several planets, 
whose orbits could be constrained by our method provided that 
the polarization fluctuations are not affected by a spectral blend 
of independent features at the same radial velocity and that the 
VLBI observations are free from a spatial blend of independent 
features along the same line of sight. The estimated mass loss 
of RLeo, ~ lO"^M0/y, is too weak to measure orbit evolu- 
tion within a reasonable amount of time (~ 3 yuas in 1000 years). 
The situation may change if the opposite effect due to tidal fric- 



tion, neglected in this estimate, is by several orders of magnitude 
more important than that of stellar mass loss. 

Without further observational evidence at hand, one can only 
speculate on whether the suggested Jovian magnetic fields are 
decisive for evaporation or survival of the planets. They sup- 
press turbulent viscosity and therefore the tidal friction induced 
by the latter (Zahn, 1977), and also lead to spin angular momen- 
tum loss due to magnetic braking. This contrasts with scenarios 
where the coupling between the stellar magnetic field and the 
planetary one lead to an increased drag area and therefore infall 
of the planet as compared to the case of an unmagnetized planet. 
The finding that our observations most likely imply a planetary 
magnetopause may turn out to be important here and motivate 
studies of planetary evolution in late-type stars that include mag- 
netic star-planet interaction (omitted by Villaver & Livio, 2007, 
but included by Nordhaus & Blackman, 2006 who predict stellar 
outflows). The case for exoplanetary magnetism proposed here 
can only be strengthened by a combination of extensive simul- 
taneous polarization monitoring and repeated spatially resolved 
observations. The experimentum crucis for an unambiguous ev- 
idence of planets will be the detection of orbital motion in the 
proper motion of maser spots showing the magnetic fluctuations, 
by contemporaneous direct imaging of the suggested wake flows 
at the same velocity as the fluctuations. Such a double evidence 
would rule out ambiguities: there may be planets without wake 
flows, if the wind density is not sufficient to produce observ- 
able features, and there are radial alignments of maser spots not 
related to wake flows. The coincidence of a persisting pseudo- 
periodic magnetic fluctuation with a radial alignment of maser 
spots strongly suggests that these phenomena are associated and 
naturally explained by a precessing Jovian magnetosphere and a 
planetary wake flow, respectively. 
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Fig. A.l. Definition of the orbital parameters. The xyz coordinate 
system is shown by solid arrows (the yz plane, i.e. the drawing 
plane, is the plane of the sky), the ^rj^ coordinate system by the 
dashed arrows (the orbit is in the ^rf plane). The rj axis is the node 
line about which the orbital plane is inclined by angle and lies 
in the xy plane. O is the longitude of the ascending node. The 
star, the planet and its wake flow are schematically indicated, as 
well as the position of the ascending and of the descending node. 



Appendix A: Determination of orbital parameters 

Here we summarize the determination of orbital parameters, 
which follows the classical method of analysing spatially and 
spectrally resolved binary systems. Without further observa- 
tional evidence at hand, we only treat the case of a circular orbit 
here. The coordinate systems applied and the orbital parame- 
ters are defined in Fig. A.l. The transformation of the radius 
vector (pointing from the star to the planet) between the orbital 
plane which defines the system ^ri( and the xyz coordinate sys- 
tem (where the x axis is the line of sight to the observer, and 
where the yz plane is the sky plane) is given by 
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where r is the distance between star and planet, is the angle 
between the orbital angular momentum vector and the direction 
pointing from the star northwards (taken positive if the positive 
^ axis is tilted in direction of the positive ^ axis), and D. - 2k jT 
is the angular velocity of the planet's orbit (with period T). The 
angles O, and QAf are the three Euler angles of the transfor- 
mation between the ^rj^ and xyz coordinate systems. They are 
related to the inclination / of the orbital angular momentum vec- 
tor with respect to the line-of-sight to the observer by 



cos / = sin cos O 



(A.2) 



At = t - t() is the time lapse between the epoch of the observa- 
tion and the reference time to when the radius vector is along the 



vector pointing from the star to the ascending node (the node line 
defines the axis about which the orbital plane is inclined with re- 
spect to the line of sight). The xyz coordinate system is station- 
ary with respect to the local standard of rest (LSR). Applying 
the same coordinate transformation to the velocity vector of the 
planet vp yields 



Vp — Vt. + Vk 



- COS COS 3) sin (Q.At) - sin C) cos (QAf) ^ 

- cos sin ct) sin (QAf) - cos <1) cos (QAf) 
+ sin sin (QAf) 



(A.3) 



where vt is the proper motion vector of the star with respect to 
the LSR (in astronomical velocity definition, i.e. negative veloc- 
ity corresponds to motion towards the observer). Here vk is the 
Keplerian speed of the planet on the assumed circular orbit, i.e. 



V'K — Q.r 



(A.4) 



where M» is the stellar mass. For the precision required here, it 
is suflicient to assume that the center of gravity of the star-planet 
system is undistinguishably close to the center of the star 

The analysis of the observations is straightforward now: We 
need to identify a predominantly radial alignment of SiO maser 
spots at the velocity where the fluctuations of circular polariza- 
tion are observed, suggesting that they are located in the wake 
flow of a planet whose magnetic dipole axis is misaligned with 
its rotation axis. The velocity at which the fluctuations occur is 
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interpreted as the line-of-sight velocity of the planet in the LSR, 
(vp)x. A sine wave fit to the measurements of 

(vp)x - (vt)x + VK-(cos0cos3)sin(QAf) - sinOcosQAf 

= (vt)x + V.sin(nAf + iA) (A.5) 

at several epochs yields the stellar LSR velocity (vt)x and the 
angular frequency Q., as well as the parameters 

y = vk Vl - cos^ O sin^ / and 
tani^ = tan <1) sec/. (A. 6) 

These parameters are used together with the positions of the 
wake flows in the VLBI maps to constrain the radius of the orbit, 
r, and therefore the Keplerian velocity v'k = 0.r, as well as the 
longitude of the ascending node, and the inclination of the or- 
bit with respect to the sky plane, from Eg. I A. 11 The stellar mass 
is then determined from lA.4] As positions, we use the head (i.e. 
the substellar tip) of the wake flow which we identify with the 
planet's position. The corresponding velocity is the maximum 
speed detected at the head of the wake flow, since the AGB wind 
is decelerated downstream the wake flow. 

The main uncertainty of this combined kinematical and as- 
trometrical method is the association between a magnetic fluc- 
tuation and a precessing planetary magnetosphere, and between 
the planet's position and the head of the wake flow. It is there- 
fore advisable to over-determine the velocity and position fit, in 
order to minimize the observational uncertainties inherent to the 
method. 



